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Abstract 



We present a new measurement of J/V' production in Pb-Pb collisions at 158 GeV/nucleon, 
from the data sample collected in year 2000 by the NA50 Collaboration, under im- 
proved experimental conditions with respect to previous years. With the target sys- 
tem placed in vacuum, the setup was better adapted to study, in particular, the most 
peripheral nuclear collisions with unprecedented accuracy. The analysis of this data 
sample shows that the (J/V') /Drell-Yan cross-sections ratio measured in the most pe- 
ripheral Pb-Pb interactions is in good agreement with the nuclear absorption pattern 
extrapolated from the studies of proton-nucleus collisions. Furthermore, this new mea- 
surement confirms our previous observation that the / Drell-Yan cross-sections 
ratio departs from the normal nuclear absorption pattern for semi-central Pb-Pb col- 
lisions and that this ratio persistently decreases up to the most central collisions. 
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1 Introduction 



The suppression of the i/ip yield in heavy ion coUisions, predicted by Matsui and 
Satz pp, is commonly considered as one of the most interesting signals of the formation 
of a deconfined state of quarks and gluons in high-energy heavy-ion collisions. The 
detection of J/?/' and ip' mesons through their leptonic decay to a pair of muons is 
particularly interesting since muons are not affected by the strong interactions at play 
in the later stages of the collision evolution, when the light hadrons are formed. 

The NA50 experiment is a high luminosity fixed target experiment at CERN, es- 
sentially dedicated to the study of dimuon production in Pb-Pb collisions at 158 GeV 
per nucleon. The analysis of the Pb-Pb data collected by NA50 in year 1995 showed j3 
Ej that the i /ifj production yield, with respect to the production of Drell-Yan dimuons, 
is ( "anomalously" ) suppressed with respect to the "normal nuclear absorption" pat- 
tern derived from measurements done in interactions induced by protons or light 
nuclei |110ini- This "integrated" result was complemented by detailed studies of the 
i/ip suppression pattern as a function of the centrality of the collision [ZIIHIE], with 
data collected in the years 1996 and 1998, which indicated that the extra suppression 
sets in for semi-central collisions and suggested that the departure from the normal 
absorption curve was setting in over a narrow range of centrality values. 

A more detailed analysis of the Pb-Pb data revealed that peripheral interactions 
could be contaminated by Pb-air interactions, especially in the multi-target config- 
uration used in the data taking period of 1996. We devoted the year 2000 data 
taking period to investigate further whether Pb-Pb peripheral collisions were really 
compatible with the results obtained from lighter collision systems, and to collect 
more statistics on several p-A systems, in order to establish a more precise normal 
absorption curve. It was also realized that, although the pattern describing their de- 
pendence as a function of the centrality of the collision remains essentially unaffected, 
the numerical values of the ratio of cross-sections B^^a{^ /ip) / (j{DY) are somewhat 
sensitive to the specific parametrization used for the Parton Distribution Functions 
(PDFs). This happens because different PDFs give slightly different shapes for the 
Drell-Yan dimuon mass distribution and, therefore, as a result of the constraint im- 
posed by the data in the high mass region, lead to a different Drell-Yan yield in the 
region 2.9 < M < 4.5 GeV/c^. Since the GRV 92 LO [lOj PDFs were used for the 
analyses of the S-U data and of the 1995 Pb-Pb data, while the MRS A (Low Q^) ^T] 
PDFs (which take into account the isospin asymmetry in the quark sea, m 7^ (i, as 
determined by the NA51 experiment were taken for later data analyses, there is 
a systematic discrepancy between the corresponding results. To overcome this small 
inconsistency, all the data referred to in this paper, both the Pb-Pb event samples 
and the lighter interacting systems used as reference, have been analyzed with a sin- 
gle set of PDFs. We have chosen the GRV 94 LO [12] PDFs, to take into account 
the quark sea asymmetry and to be consistent with the fact that the Monte Carlo 
event generators that we use are based on leading order calculations. Finally, we 
also report results normalized to the Drell-Yan yield in the muon pair mass range 
4.2 < M < 7.0 GeV/c^, which have the advantage of being essentially insensitive to 
the PDFs chosen for the analysis. 
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In summary, this paper presents the measurement of the (J/'?/') / Drell-Yan cross- 
sections ratio, as a function of coUision centrahty, from data collected in year 2000, 
using our most recent data selection and analysis procedures. Preliminary results of 
this analysis have been previously presented in Ref. jHj. 

2 Experimental setup, event reconstruction and 
data selection 

The experimental setup used by NA50 in the year 2000 included important improve- 
ments with respect to the one used in the previous Pb-Pb data taking periods of 1995, 
1996 and 1998. The main change was the introduction of a new target system under 
vacuum [1^] (see Fig.[l}, which allowed a better rejection of out-of-target interactions 
and, in particular, of Pb-air interactions. For data selection, a new target identifi- 
cation algorithm was used, based on the silicon Multiplicity Detector [TH] instead of 
the previous one, based on a system of quartz blades ^7j. Furthermore, in order 
to identify almost simultaneous multiple interactions, a new method was developed 
based on the shape analysis of the signals of the Electromagnetic Calorimeter. Except 
for the new target system, the setup is described in detail in Refs. |21IIH1- Its main 
features are recalled hereafter. 

The NA50 apparatus consists essentially of a set of beam and anti-halo counters, 
three centrality detectors, and a muon spectrometer. The intensity of the incoming 
ion beam is counted by a Beam Hodoscope (BH) made of 16 quartz slabs. Interac- 
tions occurring upstream of the target are vetoed by dedicated scintillator counters. 
Interactions in the target can be tagged by two quartz counters or by the two planes 
of silicon micro-strips of the Multiplicity Detector (MD), depicted in Fig. ^ which 
also count the charged particles produced in the angular range 1.9 < rj < 4.2. The 
other two centrality detectors are the Electromagnetic Calorimeter (EMC), which in- 
tegrates the flux of neutral transverse energy Et in the angular domain 1.1 < rj < 2.3, 
and the Zero-Degree Calorimeter (ZDC), which measures -Ezdc? ths forward energy 
carried by the beam spectator nucleons [i] > 6.3). The muon spectrometer starts with 
a hadron absorber, made of beryllium oxide, carbon and iron, followed by two sets of 
multi-wire proportional chambers, located upstream and downstream of an air-core 
toroidal magnet. The dimuon trigger is provided by 4 scintillator hodoscopes. A 
minimum bias (MB) trigger is defined by a minimal energy deposited in the ZDC 
and a "beam" trigger is provided by the BH, leading to a sample of events only con- 
strained by the detection of an identified incident Pb ion. Data in year 2000 were 
taken with the 158 GeV/nucleon beam at 1-1.4-10'^ ions/s, over 5 s bursts, and with 
a 4.0 mm thick Pb target, corresponding to 10 % of an interaction length (interaction 
probability: 9.6%). In 35 days, we collected 135 million dimuon triggers on tape 
(plus additional triggers without target or at lower beam intensities) which led to 
64 million reconstructed dimuons, of which 720 000 opposite-sign pairs in the J/ip 
mass region (2.9-3.3 GeV/c^). After data selection and background subtraction (see 
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Figure 1: The NA50 target region (beam entering from bottom right), including the 
target vacuum tank, the Muhiphcity Detector and the upstream face of the carbon 
absorber; the Electromagnetic Calorimeter (not shown) with the BeO preabsorber in 
its center is inserted between the MD and the carbon absorber during data taking. 
The inset shows part of the target vacuum tank and the upstream face of the Multi- 
plicity Detector's first plane (beam entering from left). During data taking conditions 
the target system is moved down-stream, almost touching the MD. 

below) the number of useful dimuons in the ^/ip mass region was between 100 000 
and 130 000 (depending on the particular definition of selection cuts). 

The events were processed offline through a technically improved version of the 
software, featuring a higher track reconstruction efflciency ^3 1201 with respect to the 
version used for data taking periods before year 2000. The data selection proceeded 
as follows. First of all, parasitic interactions of the incident Pb ion in the BH were re- 
jected using auxiliary scintillator counters. Then, multiple interactions were rejected 
by a temporal analysis of the signal in the EMC, allowing us to retain events where 
either one or two incident ions were detected in the BH, within a given time window, 
but only one ion interacted in the target. Some residual interaction pile-up events 
were further rejected by a diagonal band cut on the Et-Ezbc correlation, using the 
method described in Ref. [21] (Fig. |2I). 

The location of the primary interaction is determined by requiring the appropriate 
correlation between hits on the first and second planes of the MD. This method works 
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Figure 2: Correlation between transverse and forward energy, showing ±2a and ±3(T 
diagonal bands. Contours connect pixels of equal event density, ranging from 20 
(external one) to 5024 (closest to the highest density), scaled with a constant factor. 

for Et > 3 GeV, and reaches full efficiency at ~ 15 GeV (see also Ref. |IHI). It 
retains for further analysis significantly more peripheral events than the traditional 
target identification algorithm, based on quartz Cherenkov counters surrounding the 
target. 

Muons originating downstream from the target (e.g. due to collisions in the BeO 
preabsorber or in the ZDC) are rejected by a cut on p-Dxarg, where p is the momentum 
and -Dxarg is the transverse distance between the extrapolated muon track at the 
target position and the beam axis. For the analysis of this paper, which is focused 
on high mass dimuons, the level of this "muon target cut" was set at 2 % of 
probability for each muon (see Ref. |22] for details). We have verified that our results 
are stable when we increase the level of the cut up to 10 %, while of course the 
number of accepted dimuons is consequently reduced. The rather low level of the cut 
with respect to previous data taking periods is allowed by the improved experimental 
conditions in the year 2000 run, in particular the vacuum pipe extending inside the 
BeO preabsorber. We reach a quite effective rejection of the out-of-target events, 
which are especially numerous in the most peripheral centrality bin at masses around 
2.5 GeV/c2 (see also Ref. [Hj). 

Once the event selection is done, we can proceed with the definition of centrality 
classes. For this purpose we have used our three centrality detectors, with the corre- 
sponding variables: neutral transverse energy E^, forward energy Ezbc and charged 
particle multiplicity per unit of pseudo-rapidity at mid-rapidity, (dA^ch/d?7)|max, noted 
as Nch in the following for simplicity of notation j2S|- While E^ and N^h are more 
directly correlated with the energy density of the collision, Ezdc is a good estimator 
of the number of participant nucleons . The centrality classes used in this paper 
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Table 1: Centrality classes based on the transverse energy measurement. For each 
class we list the E'y range and average, together with the average and rms values of 
iVpart, h and L. 
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Table 2: Centrality classes based on the forward energy measurement. 



are listed in Tables [TJ |21 and El The bins are equidistant in the measured centrality 
estimator being used, except for the first (or first three in the case of Nch) and last 
bin. These tables include the corresponding average and rms values of the number of 
participant nucleons, iVpart, of the impact parameter, 6, and of the average length of 
nuclear matter traversed by the (pre- resonant) charmonium state, L, all of them eval- 
uated through a detailed Glauber calculation (see Ref. [21] for the Glauber formalism 
and Ref. [23] for an example of such an evaluation). For the Pb nuclear density we 
have used a 2-parameter Fermi distribution with a half-density radius of 6.624 fm and 
a diffuseness parameter of 0.549 fm (see Ref. [211 and reference Ja73 therein). The 
parameters quoted above only describe the proton distribution inside the Pb nucleus, 
and it is well known that neutrons are distributed differently. Since our reference, 
the Drell-Yan process, is isosp in-dependent, we have used a recent experimental mea- 
surement ^27j of the "neutron halo" in heavy nuclei to define a 2-parameter Fermi 
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Table 3: Centrality classes based on the measurement of the charged particle multi- 
phcity. 

distribution for neutrons, with the same half-density radius as for protons and a dif- 
fuseness parameter of 0.667 fm (for more details see Ref. [2H])- The same model was 
adopted for the U nucleus with a half-density radius of 6.8054 fm and diffuseness 
parameters of 0.605 fm for protons and 0.786 fm for neutrons. The deformation of 
the U nucleus, which interacts with random spatial orientation, was not taken into 
account. 

3 Analysis method 

Particle yields have been extracted by fitting the dimuon mass spectrum, for each 
centrality bin. First, a fit is performed to the /i^/i^ and mass spectra, to deter- 

mine the combinatorial background, mostly due to tt and K decays. The background 
in the opposite-sign dimuon sample is then parameterized according to the relation 
Nbg = 2 - R - VN++ ■ N~-, where R can be taken as unity due to the absence of 
significant charge correlations in Pb-Pb collisions. The validity of this relation is 
granted by a fiducial cut which rejects muon pairs where any of the muons would not 
have been accepted if it had the opposite charge. 

A multi-step fit to the ^~ mass spectra (Fig. El shows an example) is then 
performed to extract the four signal contributions which are relevant for masses above 
2.5 GeV/c^, namely J/^, ip' ., Drell-Yan and open charm. Finally, the acceptance 
corrections are applied, thus obtaining the cross- sect ions ratio B^^a{i /ip) / a{DY). 
This ratio of two dimuon processes has the advantage of being insensitive to absolute 
normalizations (integrated luminosity) and to most of the experimental efficiencies, 
allowing us to directly compare these values to the normal nuclear absorption pattern 
derived from lighter collision systems. Furthermore, the Drell-Yan process is known 
(see e.g. Fig. 1 in Ref. [3] and Fig. 3 in Ref. [21]) to scale linearly with the product 
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A X B of the projectile and target mass numbers, and therefore hnearly with the 
number of binary nucleon-nucleon colhsions. 

The analysis is performed in the dimuon kinematical domain < ?/cms < 1 and 
—0.5 < cos(^cs) < 0.5, where i/cms is the c.m.s. rapidity and ^cs is the polar decay 
angle of the muons in the Collins-Soper reference frame PU] . 

For each of the four signal sources, muon pairs are generated and propagated 
through the NA50 detector using the NA38/NA50 Monte Carlo simulation package 
which implements the multiple scattering and energy loss processes of muons as done 
in GEANT 3.21 (see e.g. These events are then reconstructed using the stan- 

dard NA50 reconstruction code in order to obtain the functional forms of each signal, 
needed for the fits to the measured invariant mass distributions, and their accep- 
tances, needed for the calculation of the cross-sections ratio. For J/ip and ip' genera- 
tion, we use a gaussian rapidity distribution with an r.m.s. value a = 0.6 units [3T] and 
a uniform cos(6'cs) distribution. The transverse momentum is generated according 
to Ki{Mt/T), where Ki is the modified Bessel function and T = 236 MeV [21]. The 
Drell-Yan and open charm (DD) contributions are calculated with PYTHIA 6.124 
using the GRV 94 LO ^ set of PDFs, as provided by the PDFLIB package In 
order to reproduce the measured distributions, the width of PYTHIA's gaussian 
primordial fcx distribution has been set to 0.8 GeV/c for Drell-Yan and 1.0 GeV/c 
for open charm |34j . 




Figure 3: Opposite-sign dimuon invariant mass spectrum for the centrality class 
35 < Et < 45 GeV. The results from the fit are also shown. 
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The simulated and reconstructed dimuon mass distributions are parameterized 
by empirical functional forms. Such functions are illustrated in Fig. |2l after adjust- 
ment to the measured invariant mass spectrum. A first approximation of the ^/if) 
and if)' functional forms is obtained through the Monte Carlo procedure outlined 
above, which is meant to describe the detector momentum resolution, and the mul- 
tiple scattering and energy loss in the materials (target and muon filter) crossed by 
the muons. However, the very high statistical precision of our data is incompatible 
with the imperfections of this Monte Carlo description. The calculated line shape of 
the ^/ip must be slightly adjusted in order to closely reproduce the measured J/-?/' 
peak [in|l2D|- 

The line shape of the measured ^/ip resonance has been recently studied in de- 
tail [2ni! resulting in the use of a functional form which significantly improves the fit 
of the dimuon mass distributions. The new description of the ^/ip shape leads to a 
decrease in the level of the (J/^/') / Drell-Yan cross-sections ratio by about 6% with 
respect to our preliminary results ^3], while the pattern as a function of centrality 
remains unchanged (for further details, see Refs. J3] and [SH])- 

The fit to the invariant mass distributions is made, using the maximum likelihood 
method, through a multi-step procedure. First, the combinatorial background in the 
sample is estimated through a fit to the like-sign invariant mass spectra; then a 
fit in the region 2.9 < M < 8.0 GeV/c^ is performed to obtain a first approximation 
of the J/ip, ip' and Drell-Yan contributions; then the open charm contribution is 
estimated by a fit in the mass region 1.7 < M < 2.6 GeV/c^. Finally, in order 
to minimize potential fit-induced distortions of the highly sensitive Drell-Yan and 
ip' components of the mass spectrum, the last fit is performed again in the region 
2.9 < M < 8.0 GeV/c^ (see Fig. El) with five free parameters, which are the J/'?/', 
%l>' and Drell-Yan normalizations, the ^/ip mass and the ^/ip experimental width 
parameter a^. The fitted values are in the range 97-100 MeV/c^ for cx^, and in the 
range 3.110-3.112 GeV/c^ for M^. 

The acceptances of the NA50 detector, as calculated by our Monte Carlo simula- 
tions and in the phase space window mentioned above, are given in Tabled for ^/ip 
and Drell-Yan dimuons (for the two mass regions considered in the present analysis) . 
The relative systematic error on the acceptances is estimated to be 2 %, while the 
statistical error is negligible. 



Process 


Acceptance (%) 


Drell-Yan (2.9-4.5 GeV/c^) 
Drell-Yan (4.2-7.0 GeV/c^) 


12.5 
13.8 
17.8 



Table 4: Acceptances of the experimental setup used in year 2000, for ^/ip — > /i^/i 
and Drell-Yan dimuons, in the phase space window analyzed in this study. 
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4 Results 



4.1 Summary of p-A results 

In order to accurately establish the normal nuclear absorption pattern, data were 
collected, with the NA50 apparatus, using the CERN SPS 400 or 450 GeV proton 
beam and several different nuclear targets. Results from the first two of our three 
data samples have been recently published |^ IT^ . A publication concerning the 
results from the third data sample, taken at 400 GeV beam energy, is presently 
under preparation (see Ref. for preliminary results). 

In a previous publication [21] we presented the analyses of several data sets for 
the / Drell-Yan cross-sections ratio. That study merged measurements done 

with 450 GeV protons (pp and pd from NA51; first 450 GeV p-A data samples 
from NA50), in the phase space domain —0.5 < coshes < 0.5 and —0.4 < |/cms < 0.6, 
with S-U measurements done at 200 GeV/nucleon (from NA38), in the domain 
—0.5 < cos6'cs < 0.5 and 0.0 < ?/cms < 1-0. For consistency reasons, we have revis- 
ited that study and updated its results, including now the three p-A data samples of 
NA50, to establish the normal nuclear absorption curve needed for the present work. 
Moreover, S-U results are not included anymore in this study. It is very interesting to 
check if the measured S-U (J/ip) / Drell-Yan results will become significantly different 
from the expected behaviour deduced using exclusively p-A measurements (see next 
subsection). To bring the analysis of those p-A data sets (and also of the S-U data 
set) in tune with the procedures used in the present study of the Pb-Pb data, the J/ip 
and ip' resonances, in the fits of the mass spectra, were described with the new line 
shapes [ini 1201 we consistently used the GRV 94 LO PDF sets in all calculations. 
Moreover, the Monte Carlo simulation of J/ip production in p-A collisions has been 
upgraded in order to reproduce more accurately the measured pt and y distributions 
and thus obtain more precise J/ip acceptance values for each individual data sample 
(for more details, see Ref. j2Hl)- 

The cross-sections ratios (J/V') / Drell-Yan for p-A collisions (including pp and pd 
results from NA51) are then fitted to a Glauber description of the normal nuclear 
absorption. The Glauber calculation is performed with three free parameters: the 
two independent normalizations (accounting for the different energy and kinematical 
conditions) and a common value for the J/ip absorption cross-section. The fit leads 
to (Tabs = 4.18 ± 0.35 mb, a value very similar to our previous result, 4.4 ± 0.5 mb 
(see Ref. [21] )• Energy, rapidity domain and isospin corrections are then applied to 
scale down the normalization of the absorption curve from a pp system at 450 GeV 
to a Pb-Pb system at 158 GeV. To perform this scaling, we first compare the J/ip 
absolute p-A cross-sections measured by NA50 at 450 and 400 GeV with the J/ip 
absolute p-A cross-sections measured by NA38 and NA3 j33 at 200 GeV, in the 
phase space domain —0.5 < cos(6'cs) < 0.5 and 0.0 < ?/cms < 1-0. If these absolute 
cross-sections are fitted with a common dabs, we obtain a value of 4.11 ± 0.43 mb, 
in excellent agreement with the value of 4.18 ± 0.35 mb quoted above and measured 
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using / Drell-Yan ratios at 450 and 400 GeV^. The same fit leads to a factor 

of 0.319 ± 0.025 (0.343 ± 0.027) which allows to scale down the i /ip absolute cross- 
section from the 450 (400) GeV kinematical domain to the 200 GeV kinematical 
domain. In addition, an energy correction factor of 0.737 ± 0.006 is used to scale 
down the ^/ip cross-section from 200 to 158 GeV. This correction was obtained using 
the parametrization a^{y/s) = cto(1 — M^/a/s)", with the value n = 12.8 ± 0.3 
derived from a fit to available measurements. Furthermore, the change induced by 
the different on the xp window (for a fixed ycms window: 0.0 < ?/cms < 1-0) and 
on the xp distribution (parametrized as in Ref. [35]) implies a very small correction: 
1.020 ± 0.013. The corresponding scaling factors for the Drell-Yan cross-section are 
computed at leading order using the GRV 94 LO PDF sets. The simultaneous change 
in energy and kinematical domains from 450 (400) GeV to 200 GeV amounts to 
0.504ib0.012 (0.544ib0.010), while the change in energy and xp for the much smaller 
drop between 200 and 158 GeV gives 0.7085 and 1.091, respectively. Finally, a factor 
0.969 is applied to the absorption curve corresponding to the isospin change from pp 
to Pb-Pb. 

The statistical uncertainty on the absorption curve as measured at the higher ener- 
gies is ±1 % at low centrality reaching ±4 % at high centrality. Adding in quadrature 
the errors on the rescaling factors we obtain at 158 GeV an uncertainty of ±8.3 % 
at low centrality reaching ±9.0 % at high centrality. This normal nuclear absorption 
curve is reported in the following figures along with the Pb-Pb / Drell-Yan data 

points. 

4.2 Pb-Pb results 

We now present the J/ip suppression pattern for Pb-Pb interactions, as obtained from 
the data collected in year 2000. As in the past, we consider the cross-sections ratio 
(J/ip) /Drell-Yan, where the Drell-Yan differential cross-section is integrated in the 
mass domain 2.9-4.5 GeV/c^. The results from three independent analyses, using 
the centrality variables E-p, Ezbc and A^ch, are shown in Figs. El El and El which 
also include the normal nuclear absorption curve, obtained as explained above. The 
corresponding numerical values are reported in Table El The results referred to the 
high mass Drell-Yan domain (4.2-7.0 GeV/c^), which are less sensitive to the specific 
PDF sets used in the analysis, can be obtained by scaling up those reported in Table El 
with the factor 7.96, the ratio between the Drell-Yan cross-sections integrated in the 
two mass domains. 

FigureEl(left) shows that the / Drell-Yan ratio in peripheral Pb-Pb collisions 

is perfectly consistent with the pattern of normal nuclear absorption, as deduced 
from p-A collisions alone, with the appropriate normalization to 158 GeV/nucleon 
and to Pb-Pb isospin content. The departure from the normal absorption pattern 
at ~ 35 GeV and the non saturation at high E^, already observed in previously 
published analyses, is also seen in this new sample of Pb-Pb collisions. The same 

-'^In the following, we adopt the value 4. 18 ±0.35 mb which is insensitive to most of the systematic 
uncertainties and, therefore, can be considered as our best estimate of (Tabs- 
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Figure 4: The /Drell-Yan cross-sections ratio as a function of transverse en- 

ergy for the Pb-Pb 2000 data sample (left). The normal absorption curve is presented 
together with the combined error from the Glauber fit and the rescaling procedure 
(dashed curves). The inset shows the ratio Measured/Expected, i.e. data over nor- 
mal nuclear absorption. The right panel presents the same data compared to the 
absorption curve computed using also S-U NA38 data for Uabs determination and 
curve normalization. 




(J(abs) = 4.18 mb (GRV 94 LO) 



10 15 20 25 30 35 

EzDc(TeV) 



Figure 5: The (J/'*/') /Drell-Yan cross-sections ratio as a function of forward energy 
for the Pb-Pb 2000 data sample. The inset shows the ratio Measured/Expected, i.e. 
data over normal nuclear absorption. 
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Table 5: Cross-sections ratio {J/ip) /Drell-Yan, referred to the 2.9-4.5 GeV/c^ Drell- 
Yan mass domain, for the Pb-Pb 2000 data, with three independent centrahty es- 
timators. Values for the Drell-Yan mass domain 4.2-7.0 GeV/c^ may be obtained 
as explained in the text. Classes are defined in Tables Q |21 and El Only statistical 
errors are given; systematic errors are negligible in comparison to statistical ones, as 
explained in the discussion session. 

observations can be made about the J/iJj absorption pattern as a function of the 
second centrality variable, Ezbc, as shown in Fig. El 

In Fig. m (right) we compare the Pb-Pb data to our previous determination of 
the absorption curve which made use of the NA38 S-U (J/^) /Drell-Yan ratios at 
200 GeV together with the best estimates of the p-A data (either using J/ip or 
(J/?/') / Drell-Yan results) at higher energies jSHj- This previous determination led to 
o'abs = 4.18 ± 0.35 mb. It is truly remarkable that this is exactly identical to the 
value presented above, and which was obtained through a rather different procedure, 
exclusively based on the 400 and 450 GeV proton-nucleus data. However, while the 
two determinations of (Tabs coincide, the new method leads to an absorption curve 
with a significantly larger error band (as seen in FigjH left) originating from the 
uncertainty on the energy-rapidity rescaling factors, which are no longer constrained 
by the 200 GeV S-U results. 

The analysis as a function of the third centrality variable, Nch, has the unique 
feature that the charged particle multiplicity is evaluated in the rapidity window 
where the dimuons are measured by the muon spectrometer. From the / Drell- 

Yan results displayed in Fig. El we observe that the suppression pattern is very similar 
to the one obtained with the Eq^ centrality variable. Note that the most peripheral 
data point corresponds to about half the statistics of the semi-central and central 
data points. 

A comparison of the three independent analyses is presented in Fig. [71 using the 
average number of participant nucleons, A^part, as a common centrality variable. It is 
worth recalling that -Et and the charged particle rapidity density at mid-rapidity are 
directly related to the energy density reached in the collision, through the Bjorken 
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Figure 6: The / Drell-Yan cross-sections ratio as a function of charged multiph- 

city for the Pb-Pb 2000 data sample. The inset shows the ratio Measured/Expected, 
i.e. data over normal nuclear absorption. 



formula, while the forward energy measured in the ZDC is more strongly correlated to 
the geometry of the collision, being a simple and robust estimator of A^^part- We have 
chosen iVpart for the purpose of comparing our three independent analyses because it 
is well known that, at SPS energies, both and N^-t, are linearly proportional to 
-^part, as expected in the framework of the wounded nucleon model, up to the most 
central collisions (see, for instance, Refs. j21] and j2S])- 

The horizontal error bars in Fig. [7| represent the r.m.s. values of the A^part dis- 
tribution in that bin, which depend on the experimental smearing, specific of each 
centrality variable. We see from this figure that the three (completely independent) 
centrality measurements give a very consistent picture. In particular, we do not see 
any evidence that the ^/ip suppression pattern looks different when looked as a func- 
tion of variables related to particle production {E"[ and iVch) or as a function of simple 
geometry (^zdc)- 

In order to compare our Pb-Pb results with other (lighter) collision systems, we 
present in Fig. |H1 (left) the {^/ip) / Drell-Yan cross-sections ratio obtained in this anal- 
ysis together with the results obtained in S-U (NA38) and p-A (NA50) interactions. 
The right panel of this figure shows the same data points after dividing by the func- 
tion representing the normal nuclear absorption, using for the horizontal scale the 
average length L of nuclear matter traversed by the cc state. The length L has been 
evaluated as {p{r)L) / pq, where p(r) is the nuclear density, normalized for each nu- 
cleus in order to give A nucleons upon integration, and po is the average nuclear 
density, 0.17 fm~^. The figure shows how the J/ip anomalous suppression clearly 
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Figure 7: The {J/ip) / Drcll-Yan cross-sections ratio as a function of iVpart for three 
analyses of the Pb-Pb 2000 data sample, compared to (left) and divided by (right) 
the normal nuclear absorption pattern. 
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Figure 8: The /Drell-Yan cross-sections ratio vs. L, for several collision sys- 

tems, compared to (left) and divided by (right) the normal nuclear absorption pattern. 
The measured data have been rescaled to 158 GeV/nucleon and for proton-proton 
isospin. The Et data set was used for the Pb-Pb points. 
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emerges from the reference line obtained exclusively from proton-nucleus data. We 
observe also that the S-U data points closely match the absorption curve determined 
from the p-A data, leaving very little room for anomalous absorption of the J/tp in 
S-U collisions. 
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Figure 9: The (J/ip) /Drell-Yan cross-sections ratio as a function of the energy den- 
sity, for the S-U and Pb-Pb collision systems, compared to (left) and divided by 
(right) the normal nuclear absorption pattern. The measured data have been rescaled 
to 158 GeV/nucleon and for proton-proton isospin. The Ei^ data set was used for the 
Pb-Pb points. 

Finally, we present in Fig. Elthe (J/ip) / Drell-Yan cross-sections ratio obtained in 
S-U (NA38) and Pb-Pb (NA50) interactions as a function of e, the energy density 
averaged over the whole transverse area of the collision, evaluated with the Bjorken 
formula: 

_ dEx/dr/ln^ax 

ctAt ' 

where is the overlap area in the transverse plane and r is the formation time, 
assumed to be 1 fm/c. We obtained the total transverse energy scaling up the mea- 
sured (neutral) E-^ by a factor 3. The different rapidity coverages of the NA38 and 
NA50 electromagnetic calorimeters have been taken into account, knowing that the 
dET/drj distributions depend on the collision centrality. 

This figure shows that the departure from the normal nuclear absorption curve 
sets in for energy densities around 2.5 GeV/fm^, just above the values reached in 
the most central S-U collisions. The absorption curves for S-U and Pb-Pb in the left 
panel of Fig. El are slightly different because the relation between energy density and 
L (obtained from the Glauber calculation) depends on the colliding nuclei. 
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5 Discussion 



We discuss first the statistical and systematical errors of the Pb-Pb 2000 results. 
Statistical errors in a given centrality bin, coming from the fit to the mass spectrum, 
are in the range 5.5-8%, and are, of course, dominated by the Drell-Yan statistics. 
We have evaluated the systematic errors coming from several sources (for details see 
Refs. |15| 1221 EHl): (i) variation of the mass fit starting point, (ii) variation of the 
level of the muon target cut, (iii) variation of the Et^-Ezbc correlation cut from ±3a 
to ±2a, (iv) use of a counting technique using the number of signal events in the 
mass ranges 2.9-3.3 and 4.2-7.0 GeV/c^, with subtraction of the number of Drell- 
Yan events in the first mass range. The point-to-point systematic error from the 
above sources was found to be negligible when compared to the statistical error. 

Another systematic effect comes from the choice of PDFs for the Drell-Yan pro- 
cess (influencing both the mass fit and the acceptance calculation), which affects the 
overall normalization of our results while leaving unchanged the shape of the sup- 
pression pattern as a function of centrality. We have compared our standard choice 
of PDF (GRV 94 LO) with more recent PDFs computed at leading order, namely 
GRV 98 LO t39j, MRST (central gluons) LO P| and CTEQ 5L We have studied 
the effect of changing the Drell-Yan functional form while keeping the same normal- 
ization in the 4.2-7.0 GeV/c^ mass region, as imposed by our data. The maximum 
change in the Pb-Pb (J/ip) / Drell-Yan cross- sect ions ratio is 3.5% (for comparison, 
using the sets GRV 92 LO and MRS A (Low Q^) quoted in the introduction, we 
obtain a change of 0.8% and 10%, respectively). 

We now discuss briefly how the Pb-Pb 2000 results on the / Drell-Yan cross- 

sections ratio, presented above, compare with our results from previous Pb-Pb data 
taking periods. 

We consider the 1998 data sample, collected with an experimental setup similar to 
the one used in year 2000 except for the vacuum around the target, and find excellent 
agreement, as shown in the left panel of Fig.^[ after repeating the analysis with our 
most recent procedures (see Ref. |22])- The introduction of the vacuum around the 
target region has brought a considerable increase in statistics of peripheral events. 
The weighted average of the two analyses (where data of year 1998 have been rebinned 
using the year 2000 centrality classes) is presented in the right panel of Fig. HUJ 
The resulting suppression pattern suggests two different suppression regimes, one for 
peripheral and one for central reactions. The observed transition between these two 
regimes, including the smearing due to the resolution of the e.m. calorimeter which 
amounts here to 8%, extends over a range of about 15 GeV in Et- 

6 Conclusion 

We have analyzed the most recent NA50 Pb-Pb data sample, collected in year 2000 
under improved experimental conditions, which allowed us to extend the analysis of 
the {J/ip) /Drell-Yan cross sections ratio down to very peripheral interactions using 
only dimuon triggers. We have used three different centrality estimators, namely 
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Figure 10: Comparison between 1998 and 2000 results: separate /Drell-Yan 
cross-sections ratios from 1998 and 2000 (left) and weighted average (right), as a 
function of transverse energy. 



transverse energy, forward energy and charged particle multiplicity. Our analysis 
shows that peripheral Pb-Pb interactions, with impact parameter b > 8.5 fm, exhibit 
a J/i/j production yield in agreement with the normal nuclear absorption pattern, 
with (Tabs = 4.18 mb, derived from an extensive study of p-A collisions. For smaller 
impact parameter values we observe a departure from the normal absorption curve, 
followed by a persisting decrease up to the most central Pb-Pb coUisions. 
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